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Outage Probability of the Gaussian MIMO 
Free-Space Optical Channel with PPM 

Nick Letzepis, Member, IEEE and Albert Guillen i Fabregas, Member, IEEE 

Abstract 

The free-space optical channel has the potential to facilitate inexpensive, wireless communication 
with fiber-like bandwidth under short deployment timelines. However, atmospheric effects can signifi- 
cantly degrade the reliability of a free-space optical link. In particular, atmospheric turbulence causes 
random fluctuations in the irradiance of the received laser beam, commonly referred to as scintillation. 
The scintillation process is slow compared to the large data rates typical of optical transmission. As 
such, we adopt a quasi-static block fading model and study the outage probability of the channel under 
the assumption of orthogonal pulse-position modulation. We investigate the mitigation of scintillation 
through the use of multiple lasers and multiple apertures, thereby creating a multiple-input multiple 
output (MIMO) channel. Non-ideal photodetection is also assumed such that the combined shot noise 
and thermal noise are considered as signal-independent additive Gaussian white noise. Assuming perfect 
receiver channel state information (CSI), we compute the signal-to-noise ratio exponents for the cases 
when the scintillation is lognormal, exponential and gamma-gamma distributed, which cover a wide 
range of atmospheric turbulence conditions. Furthermore, we illustrate very large gains, in some cases 
larger than 15 dB, when transmitter CSI is also available by adapting the transmitted electrical power. 

I. Introduction 

Free-space optical (FSO) communication offers an attractive alternative to the radio frequency 
(RF) channel for the purpose of transmitting data at very high rates. By utilising a high carrier 
frequency in the optical range, digital communication on the order of gigabits per second is 
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possible. In addition, FSO links are difficult to intercept, immune to interference or jamming from 
external sources, and are not subject to frequency spectrum regulations. FSO communications 
have received recent attention in applications such as satellite communications, fiber-backup, 
RF- wireless back-haul and last-mile connectivity [1]. 

The main drawback of the FSO channel is the detrimental effect the atmosphere has on a 
propagating laser beam. The atmosphere is composed of gas molecules, water vapor, pollutants, 
dust, and other chemical particulates that are trapped by Earth's gravitational field. Since the 
wavelength of a typical optical carrier is comparable to these molecule and particle sizes, the 
carrier wave is subject to various propagation effects that are uncommon to RF systems. One 
such effect is scintillation, caused by atmospheric turbulence, and refers to random fluctuations 
in the irradiance of the received optical laser beam (analogous to fading in RF systems) [2-4]. 

Recent works on the mitigation of scintillation concentrate on the use of multiple-lasers and 
multiple-apertures to create a multiple-input-multiple-output (MIMO) channel [5-13]. Many of 
these works consider scintillation as an ergodic fading process, and analyse the channel in terms 
of its ergodic capacity. However, compared to typical data rates, scintillation is a slow time- 
varying process (with a coherence time on the order of milliseconds), and it is therefore more 
appropriate to analyse the outage probability of the channel. To some extent, this has been done 
in the works of [6, 10, 12-14]. In [6, 13] the outage probability of the MIMO FSO channel is 
analysed under the assumption of ideal photodetection (i.e. a Poisson counting process) with no 
bandwidth constraints. Wilson et al. [10] also assume perfect photodetection, but with the further 
constraint of pulse-position modulation (PPM). Lee and Chan [12], study the outage probability 
under the assumption of on-off keying (OOK) transmission and non-ideal photodetection, i.e. 
the combined shot noise and thermal noise process is modeled as zero mean signal independent 
additive white Gaussian noise (AWGN). Farid and Hranilovic [14] extend this analysis to include 
the effects of pointing errors. 

In this paper we study the outage probability of the MIMO FSO channel under the assumptions 
of PPM, non-ideal photodetection, and equal gain combining (EGC) at the receiver. In particular, 
we model the channel as a quasi- static block fading channel whereby communication takes 
place over a finite number of blocks and each block of transmitted symbols experiences an 
independent identically distributed (i.i.d.) fading realisation [15, 16]. We consider two types of 
CSI knowledge. First we assume perfect CSI is available only at the receiver (CSIR case), and the 
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transmitter knows only the channel statistics. Then we consider the case when perfect CSI is also 
known at the transmitter (CSIT case). 1 Under this framework we study a number of scintillation 
distributions: lognormal, modelling weak turbulence; exponential, modelling strong turbulence; 
and gamma-gamma [17], which models a wide range of turbulence conditions. For the CSIR 
case, we derive signal-to-noise ratio (SNR) exponents and show that they are the product of: 
a channel related parameter, dependent on the scintillation distribution; the number of lasers 
times the number of apertures, reflecting the spatial diversity; and the Singleton bound [18- 
20], reflecting the block diversity. For the CSIT case, the transmitter finds the optimal power 
allocation that minimises the outage probability [21]. Using results from [22], we derive the 
optimal power allocation subject to short- and long-term power constraints. We show that very 
large power savings are possible compared to the CSIR case. Interestingly, under a long-term 
power constraint, we show that delay-limited capacity [23] is zero for exponential and (in some 
cases) gamma-gamma scintillation, unless one codes over multiple blocks, and/or uses multiple 
lasers and apertures. 

The paper is organised as follows. In Section II, we define the channel model and assumptions. 
In Section III we review the lognormal, exponential and gamma-gamma models. Section IV 
defines the outage probability and presents results on the minimum- mean squared error (MMSE). 
Then in Sections V and VI we present the main results of our asymptotic outage probability 
analysis for the CSIR and CSIT cases, respectively. Concluding remarks are then given in 
Section VII. Proofs of the various results can be found in the Appendices. 

II. System Model 

We consider an M x N MIMO FSO system with M transmit lasers an N aperture receiver as 
shown in Fig. 1 . Information data is first encoded by a binary code of rate R c . The encoded stream 
is modulated according to a Q-dsy PPM scheme, resulting in rate R = R c log 2 Q (bits/channel 
use). Repetition transmission is employed such that the same PPM signal is transmitted in perfect 
synchronism by each of the M lasers through an atmospheric turbulent channel and collected 
by N receive apertures. We assume the distance between the individual lasers and apertures is 

'Given the slow time-varying scintillation process, CSI can be estimated at the receiver and fed back to the transmitter via a 
dedicated feedback link. 
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sufficient so that spatial correlation is negligible. At each aperture, the received optical signal 
is converted to an electrical signal via photodetection. Non-ideal photodetection is assumed 
such that the combined shot noise and thermal noise processes can be modeled as zero mean, 
signal independent AWGN (an assumption commonly used in the literature, see e.g. [3-5, 12, 



In FSO communications, channel variations are typically much slower than the signaling 
period. As such, we model the channel as a non-ergodic block-fading channel, for which a given 
codeword of length BL undergoes only a finite number B of scintillation realisations [15, 16]. 
The received signal at aperture n, n = 1, . . . , N can be written as 



for b = 1, . . . , B, £ = 1, . . . , L, where y b [£\, z b [£] e M Q are the received and noise signals at 
block b, time instant £ and aperture n, x b [£], e M Q is the transmitted signal at block b and time 
instant £, and h b n ' n denotes the scintillation fading coefficient between laser m and aperture n. 
Each transmitted symbol is drawn from a PPM alphabet, Xb[£] e X ppm = {e±, . . . , cq}, where 
e q is the canonical basis vector, i.e., it has all zeros except for a one in position q, the time 
slot where the pulse is transmitted. The noise samples of z b [£] are independent realisations of a 
random variable Z ~ jV(0, 1), and p b denotes the received electrical power of block b at each 
aperture in the absence of scintillation. The fading coefficients /i™' n are independent realisations 
of a random variable H with probability density function (pdf) ffj(h). 

At the receiver, we assume equal gain combining (EGC) is employed, such that the entire 
system is equivalent to a single-input single-output (SISO) channel, i.e. 



where z b [£] = ^E n =i z fcM ~ A/"(0, 1), and h b , a realisation of the random variable H, is 
defined as the normalised combined fading coefficient, i.e. 

M N 



14, 24-29]). 




(1) 




n=l 



(2) 




(3) 



m=l n=l 
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where c = l/(M[H]y/l + a]/(MN)) is a constant to ensure E[H 2 } = l. 2 Thus, the total 
instantaneous received electrical power at block b is p b — M 2 Np b /c, and the total average 
received SNR is snr = E[/if,p&] = E[p&]. 

For the CSIR case, we assume the electrical power is distributed uniformly over the blocks, i.e., 
p b — p — snr for b — 1, . . . , B. Otherwise, for the CSIT case, we will allocate electrical power 
in order to improve performance. In particular, we will consider the following two electrical 
power constraints 

1 B 

Short-term: — ^p b < P (4) 



B 

6=1 



Long-term: E 



1 A 



B 

6=1 



< P- (5) 



Throughout the paper, we will devote special attention to the case of B = 1, i.e., the channel 
does not vary within a codeword. This scenario is relevant for FSO, since, due to the large 
data-rates, one is able to transmit millions of bits over the same channel realisation. We will 
see that most results admit very simple forms, and some cases, even closed form. This analysis 
allows for a system characterisation where the expressions highlight the roles of the key design 
parameters. 

III. Scintillation Distributions 
The scintillation pdf, fjj(h), is parameterised by the scintillation index (SI), 

2 A Var(^) 

Oj = = — , (o) 

1 {E[H]) 2 

and can be considered as a measure of the strength of the optical turbulence under weak 
turbulence conditions [17,30]. 

The distribution of the irradiance fluctuations is dependent on the strength of the optical 
turbulence. For the weak turbulence regime, the fluctuations are generally considered to be 

2 For optical channels with ideal photodetection, the normalisation K[H] = 1 is commonly used to keep optical power constant. 
We assume non-ideal photodetection and work entirely in the electrical domain. Hence, we chose the normalisation E[ff 2 ] = 1, 
used commonly in RF fading channels. However, since we consider only the asymptotic behaviour of the outage probability, 
the specific normalisation is irrelevant and does not affect our results. 
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lognormal distributed, and for very strong turbulence, exponential distributed [2,31]. For mod- 
erate turbulence, the distribution of the fluctuations is not well understood, and a number of 
distributions have been proposed, such as the lognormal-Rice distribution [4, 17, 32-34] (also 
known as the Beckmann distribution [35]) and K-distribution [32]. In [17], Al-Habash et al. 
proposed a gamma-gamma distribution as a general model for all levels of atmospheric tur- 
bulence. Moreover, recent work in [34] has shown that the gamma-gamma model is in close 
agreement with experimental measurements under moderate-to-strong turbulence conditions. In 
this paper we focus on lognormal, exponential, and gamma-gamma distributed scintillation, which 
are described as follows. 

For lognormal distributed scintillation, 

f±(h) = exp (-(log h - rf/(2a 2 )) , (7) 

hoy In 

where /i and a are related to the SI via /i = — log(l + <rf) and a 2 = log(l + erf). 
For exponential distributed scintillation 

f c ^(h) = \exp(-\h) (8) 

which corresponds to the super-saturated turbulence regime, where aj — 1. 

The gamma-gamma distribution arises from the product of two independent Gamma distributed 
random variables and has the pdf [17], 

f ^ k) = fP!^" 1 K ^( 2 V^) 5 (9) 
where K v (x) denotes the modified Bessel function of the second kind [36, Ch. 10]. The param- 
eters a and (3 are related with the scintillation index via a] = a^ 1 + + (a/?) -1 . 

IV. Information Theoretic Preliminaries 

The channel described by (2) under the quasi-static assumption is not information stable [37] 
and therefore, the channel capacity in the strict Shannon sense is zero. It can be shown that the 
codeword error probability of any coding scheme is lower bounded by the information outage 
probability [15, 16], 

P out (snr, R) = Pr(/(p, h) < R), (10) 
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where R is the transmission rate and I(p, h) is the instantaneous input-output mutual information 
for a given power allocation p = (pi, . . . ,pb), and vector channel realisation h = (hi, . . . , Kb). 
The instantaneous mutual information can be expressed as [38] 

1 3 

I(p,h) = -J2l aWgn (Pi>hl), (ID 

6=1 

where J awgn (p) is the input-output mutual information of an AWGN channel with SNR p. For 
PPM [24] 



/ awgn (p)=log 2 Q-E 



Q 

log 2 ( l + J2e- p+ ^ z "- Zl) 

q=2 



(12) 



where Z q ~ Af(0, 1) for q = 1, . . . , Q. 

For the CSIT case we will use the recently discovered relationship between mutual information 
and the MMSE [39]. This relationship states that 3 

dp [P) log(2) {U) 
where mmse(p) is the MMSE in estimating the input from the output of a Gaussian channel as 
a function of the SNR p. For PPM, we have the following result 

Theorem 4.1: The MMSE for PPM on the AWGN channel with SNR p is 



mmse(p) = 1 — E 



(14) 



where Z { ~ Af(0, 1) for i = 1, . . . , Q. 

Proof: See Appendix I. ■ 
Note that both (12) and (14) can be evaluated using standard Monte-Carlo methods. 

V. Outage Probability Analysis with CSIR 

For the CSIR case, we employ uniform power allocation, i.e. p± — . . . — p B — snr. For 
codewords transmitted over B blocks, obtaining a closed form analytic expression for the outage 
probability is intractable. Even for B — 1, in some cases, for example the lognormal and gamma- 
gamma distributions, determining the exact distribution of H can be a difficult task. Instead, as 

3 The log(2) term arises because we have defined J awgn (p) in bits/channel usage. 
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we shall see, obtaining the asymptotic behaviour of the outage probability is substantially simpler. 
Towards this end, and following the footsteps of [20,40], we derive the SNR exponent. 

Theorem 5.1: The outage SNR exponents for a MIMO FSO communications system modeled 
by (2) are given as follows: 

) 2 = siogaT^) (1 + [B (1 - ^ )J) (15) 

MN 

^Ssn^^-a+L^a-^j), (i6) 

MN 

4 g ogsnr) = ~2~ min(a,/3) (1 + [B (1 - R c )\) , (17) 

for lognormal, exponential, and gamma-gamma cases respectively, where R c = R/\og 2 (Q) is 
the rate of the binary code and 

, a logP out (snr, J R) , 

«(logsnr)fe = - lim T, rr k = 1, 2. (18) 

v ' snr^oo (logSnr) ft 

Proof: See Appendix EL ■ 

Proposition 5.1: The outage SNR exponents given in Theorem 5.1, are achievable by random 
coding over PPM constellations whenever B (1 — R c ) is not an integer. 

Proof: The proof follows from the proof of Theorem 5.1 and the proof of [20, Th. 1]. ■ 

The above proposition implies that the outage exponents given in Theorem 5.1 are the optimal 
SNR exponents over the channel, i.e. the outage probability is a lower bound to the error 
probability of any coding scheme, its corresponding exponents (given in Theorem 5.1) are an 
upper bound to the exponent of coding schemes. From Proposition 5.1, we can achieve the 
outage exponents with a particular coding scheme (random coding, in this case), and therefore, 
the exponents in Theorem 5.1 are optimal. 

From (15)-(17) we immediately see the benefits of spatial and block diversity on the system. In 
particular, each exponent is proportional to: the number of lasers times the number of apertures, 
reflecting the spatial diversity; a channel related parameter that is dependent on the scintillation 
distribution; and the Singleton bound, which is the optimal rate-diversity tradeoff for Rayleigh- 
faded block fading channels [18-20]. 

Comparing the channel related parameters in (15)-(17) the effects of the scintillation distribu- 
tion on the outage probability are directly visible. For the lognormal case, the channel related 
parameter is 81og(l + erf) and hence is directly linked to the SI. Moreover, for small erf < 1, 
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8 log(l + aj ) w 8(7j and the SNR exponent is inversely proportional to the SI. For the exponential 
case, the channel related parameter is a constant 1/2 as expected, since the SI is constant. 
For the gamma-gamma case the channel related parameter is min(ai,/3)/2, which highlights an 
interesting connection between the outage probability and recent results in the theory of optical 
scintillation. For gamma-gamma distributed scintillation, the fading coefficient results from the 
product of two independent random variables, i.e. H = XY , where X and Y model fluctuations 
due to large scale and small scale cells. Large scale cells cause refractive effects that mainly 
distort the wave front of the propagating beam, and tend to steer the beam in a slightly different 
direction (i.e. beam wander). Small scale cells cause scattering by diffraction and therefore distort 
the amplitude of the wave through beam spreading and irradiance fluctuations [4, p. 160]. The 
parameters a, (3 are related to the large and small scale fluctuation variances via a = a x 2 and 
P = ay 2 - For a plane wave (neglecting inner/outer scale effects) a Y > a 2 x , and as the strength of 
the optical turbulence increases, the small scale fluctuations dominate and a\ — > 1 [4, p. 336]. 
This implies that the SNR exponent is exclusively dependent on the small scale fluctuations. 
Moreover, in the strong fluctuation regime, a\ — > 1, the gamma-gamma distribution reduces to 
a K-distribution [4, p. 368], and the system has the same SNR exponent as the exponential case 
typically used to model very strong fluctuation regimes. 

In comparing the lognormal exponent with the other cases, we observe a striking difference. 
For the lognormal case (15) implies the outage probability is dominated by a (log(snr)) 2 term, 
whereas for exponential and gamma-gamma scintillation it is dominated by a log(snr) term. Thus 
the outage probability decays much more rapidly with SNR for the lognormal case than it does 
for the exponential or gamma-gamma cases. Furthermore, for the lognormal case, the slope of 
the outage probability curve, when plotted on a log-log scale, will not converge to a constant 
value. In fact, a constant slope curve will only be observed when plotting the outage probability 
on a log-(log) 2 scale. 

In deriving (15) (see Appendix II-A) we do not rely on the lognormal approximation 4 , which 
has been used on a number occasions in the analysis of FSO MIMO channels, e.g. [5, 12, 
29]. Under this approximation, H is lognormal distributed (7) with parameters /x = — log(l + 

4 This refers to approximating the distribution of the sum of lognormal distributed random variables as lognormal [41-44]. 
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aj/(MN)) and a 2 = — //, and we obtain the approximated exponent 

d(logsnrP ~ — (1 +[B(l~ Rc)\ ) • (19) 

Comparing (15) and (19) we see that although the lognormal approximation also exhibits a 
(log(snr)) 2 term, it has a different slope than the true SNR exponent. The difference is due to 
the approximated and true pdfs having different behaviours in the limit as h — > 0. However, for 
very small aj < 1, using log(l + x) x (for x < 1) in (15) and (19) we see that they are 
approximately equal. 

For the special case of single block transmission, B = 1, it is straightforward to express 
the outage probability in terms of the cumulative distribution function (cdf) of the scintillation 
random variable, i.e. 

( /snr awgn \ 

P out (snr, R)=F H U —L \ (20 ) 

where F H (h) denotes the cdf of H, and snr awgn = 7 awgn '^(R) denotes the SNR value at which 
the mutual information is equal to R. Table I reports these values for Q = 2,4,8,16 and 
R — R c log 2 Q, with R c — \, |, |. Therefore, for B = 1, we can compute the outage probability 
analytically when the distribution of H is available, i.e., in the exponential case for M, N > 1 or 
in the lognormal and gamma-gamma cases for M, N = 1 . In the case of exponential scintillation 
we have that 



sn| ,awgn\ j 

P out (snr, J R) = f ( MN, [MN(l + MN)^-j ], (21) 



where T(a,x) = J* t a ~ x exp(— t) dt denotes the regularised (lower) incomplete gamma 
function [36, p. 260]. For the lognormal and gamma-gamma scintillation with MN > 1, we 
must resort to numerical methods. This involved applying the fast Fourier transform (FFT) to 
ffj to numerically compute its characteristic function, taking it to the MiV th power, and then 
applying the inverse FFT to obtain f H . This method yields very accurate numerical computations 
of the outage probability in only a few seconds. 

Outage probability curves for the B = 1 case are shown on the left in Fig. 2. For the lognormal 
case, we see that the curves do not have constant slope for large SNR, while, for the exponential 
and gamma-gamma cases, a constant slope is clearly visible. We also see the benefits of MIMO, 
particularly in the exponential and gamma-gamma cases, where the SNR exponent has increased 
from 1/2 and 1 to 2 and 4 respectively. 
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VI. Outage Probability Analysis with CSIT 



In this section we consider the case where the transmitter and receiver both have perfect CSI 
knowledge. In this case, the transmitter determines the optimal power allocation that minimises 
the outage probability for a fixed rate, subject to a power constraint [21]. The results of this 
section are based on the application of results from [22] to PPM and the scintillation distributions 
of interest. Using these results we uncover new insight as to how key design parameters influence 
the performance of the system. Moreover, we show that large power savings are possible 
compared to the CSIR case. 

For the short-term power constraint given by (4), the optimal power allocation is given by 
mercury- waterfilling at each channel realisation [22,45], 



for b = 1, . . . , B where mmse -1 ^) is the inverse-MMSE function and rj is chosen to satisfy the 
power constraint. 5 From [22, Prop. 1] it is apparent that the SNR exponent for the CSIT case 
under short-term power constraints is the same as the CSIR case. 

For the long-term power constraint given by (5) the optimal power allocation is [22] 




(22) 




(23) 



where 




(24) 



and s is a threshold such that s = oo if lim^oo E^( s ) 



[i Ef=i Pb < P, and 




(25) 



otherwise, s is chosen such that P = -| J2b=i P&J • m (24), f] 

rate constraint 



is now chosen to satisfy the 




(26) 



5 Note that in [22,45], the minimum in (22) is between 1 and ^. For QPPM, mmse(O) = (see (14)). Hence we must 
replace 1 with 
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From [22], the long-term SNR exponent is given by 

^ <! , (27) 

°° > 1 

where G^ ogsnr ) is the short-term SNR exponent, i.e., the SNR exponents (15)-(17). Note that 
^(logsnr) = 00 implies the outage probability curve is vertical, i.e. the power allocation scheme (23) 
is able to maintain constant instantaneous mutual information (11). The maximum achievable 
rate at which this occurs is defined as the delay-limited capacity [23]. From (27) and (15)-(17), 
we therefore have the following corollary. 

Corollary 6.1: The delay-limited capacity of the channel described by (2) with CSIT subject 
to long-term power constraint (5) is zero whenever 

{2 (1 + [B (1 - R c )\ y 1 exponential 
(28) 
— (1 + [B(1- i? c )J) _1 gamma-gamma 

For lognormal scintillation, delay-limited capacity is always nonzero. 

Corollary 6.1 outlines fundamental design criteria for nonzero delay-limited capacity in FSO 
communications. Single block transmission (B = 1) is of particular importance given the slow 
time-vary nature of scintillation. From (28), to obtain nonzero delay-limited capacity with B = 1, 
one requires MN > 2 and MN > 2/min(o;, 0) for exponential and gamma-gamma cases 
respectively. Note that typically, a, (3 > 1. Thus a 3 x 1, 1 x 3 or 2 x 2 MIMO system is 
sufficient for most cases of interest. 

In addition, for the special case B = 1, the solution (24) can be determined explicitly since 

r] = (/ i 2 mmse(/ awsn - 1 (i?))) _1 = (/i 2 mmse(snr awgn )) _1 . (29) 

Therefore, 

awgn 

P° pt = (30) 

Intuitively, (30) implies that for single block transmission, whenever snr a ^ gn /h 2 < s, one simply 
transmits at the minimum power necessary so that the received instantaneous SNR is equal to 
the SNR threshold (snr^ 11 ) of the code. Otherwise, transmission is turned off. Thus an outage 



occurs whenever h < \ — B — and hence 



awgn 



Po,u(snr.ro= /'//h/^r^l (3D 
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where 7 1 (snr) is the solution to the equation 7(5) = snr, i.e., 

7 (a) ± snr^ n ^ ^ dh, (32) 

A /snr awgn 

where z/ = y — | — . Moreover, the snr at which P out (R, snr) — > is precisely lim^oo 7(5). 
In other words, the minimum long-term average SNR required to maintain a constant mutual 
information of R bits per channel use, denoted by snr, is 

snr77 gn = snr* wgn f°° dh. (33) 

Jo h 

Hence, recalling that snr^ 11 = J awgn, ~ 1 (i?), the delay-limited capacity (under the constraint 
of PPM) is 6 

c ^ - *~ ■ (34) 

In the cases where the distribution of H is known in closed form, (32) can be solved explicitly, 
hence yielding the exact expressions for outage probability (31) and delay-limited capacity (34). 
For lognormal distributed scintillation with B — M — N — 1, we have that 

V(.) = i-rj-d + erfc f 3 1 o g (l +J ?) + |.o g sn rr °-l logs \ 



and 



C(snr) = J awgn ( ) , (36) 



:i+^ 2 ) 

where we have explicitly solved the integrals in (32) and (34) respectively 
For the exponential case with B = 1, we obtain, 



-aw gn MN(1 + MN) - / o /„ ^ ..^srn^ 



awgn 



and 



7 -( S ) = snr r n (MiV _ V 1)(Miv : 2) r I 1 / -V - 2. y 3 / ,V( 1 1 AA') I . (37) 

/awgn / (MN-l)(MN-2) \ MN>2 

Cr(sur)={ y MN{1+MN) ) (38) 

otherwise. 



6 Note that a similar expression was derived in [23]. 
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For the gamma-gamma case with B — M — N — 1, 7 gg (s) can be expressed in terms of 
hypergeometric functions, which are omitted for space reasons. The delay-limited capacity, 
however, reduces to a simpler expression 7 



Fig. 2 (right) compares the outage probability for the B = 1 CSIT case (with long-term power 
constraints) for each of the scintillation distributions. For MN = 1 we see that the outage curve 
is vertical only for the lognormal case, since Cd = for the exponential and gamma-gamma 
cases. In these cases one must code over multiple blocks for C d > 0, i.e. from Corollary 6.1, 
B > 6 and B > 4 for the exponential and gamma-gamma cases respectively (with R c = 1/2). 
Comparing the CSIR and CSIT cases in Fig. 2 we can see that very large power savings are 
possible when CSI is known at the transmitter. These savings are further illustrated in Table II, 
which compares the SNR required to achieve P out < 1CT 5 (denoted by snr*) for the CSIR case, 
and the long-term average SNR required for P out — > in the CSIT case (denoted by snr, which 
is given by (33)). Note that in the CSIT case, the values of snr given in the parentheses' is the 
minimum SNR required to achieve P out < 1CT 5 , since Cd = for these cases (i.e. snr = oo). 
From Table II we see that the power saving is at least around 15 dB, and in some cases as high 
as 50 dB. We also see the combined benefits of MIMO and power control, e.g. at MN = 4, the 
system is only 3.7 dB (lognormal) to 5.2 dB (exponential) from the capacity of nonfading PPM 
channel (snr;| n = 3.18 dB). 



In this paper we have analysed the outage probability of the MIMO Gaussian FSO channel 
under the assumption of PPM and non-ideal photodetection, for lognormal, exponential and 
gamma-gamma distributed scintillation. When CSI is known only at the receiver, we have shown 
that the SNR exponent is proportional to the number lasers and apertures, times a channel related 
parameter (dependent on the scintillation distribution), times the Singleton bound, even in the 
cases where a closed form expression of the equivalent SISO channel distribution is not available 

7 Note that since we assume the normalisation E[H 2 ] = 1, then / °° S - M ^fi- dh = 4j- J °° - 3 Jfi- du, where c = 1 / yj 1 + a'j 
and ff(h) is denned as in (9) such that E[H] = 1. 




otherwise. 



(39) 



VII. Conclusion 
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in closed-form. When the scintillation is lognormal distributed, we have shown that the outage 
probability is dominated by a (log(snr)) 2 term, whereas for the exponential and gamma-gamma 
cases it is dominated by a log(snr) term. When CSI is also known at the transmitter, we applied 
the power control techniques of [22] to PPM to show that very significant power savings are 
possible. 

Appendix I 
Proof of Theorem 4.1 

Suppose PPM symbols are transmitted over an AWGN channel, the non-fading equivalent 
of (2). The received noisy symbols are given by y = ^fpx + z, where x E X ppm (we have 
dropped the time index £ for brevity of notation). 

Using Bayes' rule [46], the MMSE estimate is 

* = E[.|y] = f; _y^f»> ■ (40) 
q =i 2X1 exp{y/PVk) 

From (40) the ith element of a; is 

Using the orthogonality principle [47] mmse(p) = E[||x-i|| 2 ] = E[||a:|| 2 ] - E[||:e|| 2 ]. Since 
||e g || 2 = 1 for all q — 1, . . . , Q, then E[||:e|| 2 ] = 1. Due to the symmetry of QPPM we need 
only consider the case when x = e\ was transmitted. Hence, 

mmse(p) = 1 - (E[£ 2 ] + (Q - 1)E[£ 2 ]) . (42) 

Now y 1 = yfp + z\ and yi — Zi for % — 2, . . . , Q, where z q is a realisation of a random 
variable Z q ~ jV(0, 1) for q — 1, . . . , Q. Hence, substituting these values in (41) and taking the 
expectation (42) yields the result given the theorem. 

Appendix II 
Proof of Theorem 5.1 



We begin by defining a normalised (with respect to SNR) fading coefficient, (™' n = 
which has a pdf 



21og/C 
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Since we are only concerned with the asymptotic outage behaviour, the scaling of the coefficients 
is irrelevant, and to simplify our analysis we assume E[H 2 ] = 1. Hence the instantaneous SNR 
for block b is given by 

( M N \ 2 

— £5>r*(i-C")j (44) 
m=l n=l / 

for b = 1, . . . , B. Therefore, 

{0 if all C' n > 1 

log 2 Q at least one (™' n < 1 

= log 2 Q(l-i{C 6 ^l}) 

where Cb = (C&'\ • • • , C^). *{•} denotes the indicator function, 1 = (1, . . . , 1) is a 1 x MN 
vector of l's, and the notation a y b for vectors a,b £ R k means that Oj > bi for % — 1, . . . , k. 
From the definition of outage probability (10), we have 

P out (snr, R) = Pr(J h (snr) < R) = [ f(C)dC (45) 

J A 



where C — (Ci> • • • > Cb) is a 1 x BMN vector of normalised fading coefficients, /(£) denotes 
their joint pdf, and 

A = jc G R BMJV = E 11 ^ 1} > 5 (1 - i?c)| (46) 
is the asymptotic outage set. We now compute the asymptotic behaviour of the outage probability, 



i.e. 



- lim logP out (snr,i?) = - lim log / /(C)dC (47) 

A. Lognormal case 

From (7) and (43) we obtain the joint pdf, 

/(0=exp(-^£££(O»)A, (48) 

\ 6=1 m =l n=l / 

where we have ignored terms of order less than (logsnr) 2 in the exponent and constant terms inde- 
pendent of C in front of the exponential. Combining (47), (48), and using Varadhan's lemma [48], 



- lim logP out (snr 



, , 2 (BMN ^ 

'*> = fl0 S^ f EEE(C") 2 

I, fe=l m =l rt=l J 
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The above infimum occurs when any k of the C, b vectors are such that C, b y 1 and the other 
B — k vectors are zero, where k is a unique integer satisfying 

k < B(l -R c ) < k + 1. (49) 



Hence, it follows that k, — 1 + [B (1 — i? c )J and thus, 

(logsnr^ 2 

"8a 



lim log P out (snr, R) = 1 ; MiV (1 + |£ (1 - P C )J ) • (50) 



Dividing both sides of (50) by (logsnr) 2 the SNR exponent (15) is obtained. 

B. Exponential case 

From (8) and (43) we obtain the joint pdf, 

( MM B M N \ 

/(C) = exp - log snr^ ^ E E C" l (51) 

\ fe=l m=l n=l / 

where we have ignored exponential terms in the exponent and constant terms independent of £ 
in front of the exponential. 

Following the same steps as the lognormal case i.e. the defining the same asymptotic outage 
set and application of Varadhan's lemma [48], the SNR exponent (16) is obtained. 



C. Gamma-gamma case 
Let us first assume a > f3. From (9) and (43) we obtain the joint pdf, 

/cHC) = exp (-^C logsnr), C>0 (52) 

for large snr, where we have used the approximation K v (x) ps \T(i>)(\x)~ v for small x and 
v > [36, p. 375]. The extra condition, ( > 0, is required to ensure the argument of the 
Bessel function approaches zero as snr — > oo to satisfy the requirements of the aforementioned 
approximation. For the case (3 > a we need only swap a and (3 in (52). Hence we have the 
joint pdf 

/(C)=exp(- mi "'°-f° gSnr X: t t <r), <^0. (53) 

\ 6=1 m=l n=l / 

Now, following the same steps as in the lognormal and exponential cases, with the additional 
constraint £ b y 0, the SNR exponent (17) is obtained 
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Fig. 1. Block diagram of an M x N MIMO FSO system. 
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snr (dB) snr (dB) 



Fig. 2. Outage probability curves for the CSIR (left) and CSIT (right) cases with aj = 1, B = 1, Q = 2, R c = 1/2, 
snr^J| n = 3.18 dB: lognormal (solid); exponential (dashed); and, gamma-gamma distributed scintillation (dot-dashed), a = 2, 

13 = 3. 
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TABLE I 

Minimum signal-to-noise ratio snr^ 8 " (in decibels) for reliable communication for target rate 

R = R c log 2 Q. 



Q 


R c — \ 


R c —\ 


Rc — | 


2 


-0.7992 


3.1821 


6.4109 


4 


0.2169 


4.0598 


7.0773 


8 


1.1579 


4.8382 


7.7222 


16 


1.9881 


5.5401 


8.3107 
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TABLE II 

Comparison of CSIR and CSIT cases with B = 1, R = 1/2, Q = 2 a? = 1, a = 2, (3 = 3. Both snr* and snr are 

MEASURED IN DECIBELS. 





lognormal 


exponential 


gamma- 


•gamma 


MN 


snr* 


snr 


snr* 


snr 


snr* 


sTTF 


1 


40.1 


15.2 


106.2 


(56.2) 


65.6 


(24.5) 


2 


29.2 


9.9 


57.9 


(17.8) 


40.7 


12.2 


3 


24.4 


7.9 


42.0 


11.0 


31.7 


9.0 


4 


21.5 


6.9 


34.1 


8.4 


26.9 


7.5 
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